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Abstract 

A wide class of phase space distributions of a single mode radiation 
field is shown to be directly accessible to measurement by linear symmetric 
three-port optical couplers. 

In a classical description of the harmonic oscillator the measurement of some 
variable can be performed, at least in principle, with arbitrary precision. In 
particular, a pair of conjugated variables can jointly be measured, thus provid- 
ing a unique correspondence between states of the oscillator and points in the 
phase space of complex amplitudes. This is no longer true in quantum mechan- 
ics, to which classical phase space description cannot be transferred [|l]. Two 
main features have to be considered: i) the uncertainty principle which prevents 
precise identification of points in phase space and ii) the commutation relations 
which lead to ordering rules for representing physical variables in terms of oper- 
ators Nevertheless, a phase space description of the harmonic oscillator has 
been very fruitful in quantum optics and different distribution functions have 
been introduced to describe various dynamical processes. Some examples are 
laser theory, the study of the phase properties of light || and the quantum state 
measurement problem j| |E| . 

An important class of quantum phase space distributions are the so-called 
generalized Wigner functions 

W s (a,a) = J ^ Xs (A,A)e^~ Aa , (1) 

where 

Xs(\,X)=Tr{pe Xaf -~ Xa+ i s ^ 2 } ( 2 ) 

is the s-ordered characteristic function of the field p and is related to different 
ordering of the boson operators a, since the statistical average (a k a l ) s - 
employing the quasi-probability distribution W s (a,a) - provides the quantum 
expectation value of the s-ordered operator product {a^ k a l } s . 
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For s = 1 we get Glauber's 'P-function 0, for s = 0, the distribution originally 
proposed by Wigner jlj and for s = — 1 the Husimi Q- function Q which are 
related to normal, symmetric and antinormal ordering of the boson operators 
respectively. Such generalized Wigner functions can be sampled experimentally 
only for parameter s < — 1, since for greater values of s negative or singular 
probabilities occur for some quantum states. 

A more general class of phase distributions has been introduced in quantum 
optics ]|]. These are defined as a convolution of the Wigner function of the 
signal Wq\s with that of the probe Wq\p and serve as a description of the signal 
in terms of the probe: 

K S p(a,a)= [ ^ W ols (a + (3,a + P) W \p{P,P) . (3) 
After some calculation this can be shown to be equivalent to 

K SP {a, a) = ilk {psD(a)ppD\a)} , (4) 

which represents K$p in terms of the overlap between the signal field ps and the 
a-displaced probe field pp; D(ct) = exp{aa^ — aa} is the displacement operator. 
Phase space density as in Eq. (||) have been introduced to account for the effect 
of the measuring apparatus in a joint measurement of conjugated variables ||. 
More recently, they also have been used in entropic descriptions of quantum 
mechanical states Q. 

In this paper we consider triple coupler homodyning of a single mode radia- 
tion field, that is, the use of a linear three-port coupler to couple the signal beam 
with a strong local oscillator and also with a second probe field. Upon using 
vacuum as the probe mode input this leads to the measurement of Husimi's Q- 
function, a result which is a by-product of the present approach and is somewhat 
interesting since the measurement of the Q-function is commonly associated with 
an 4 x 4 - port device || (eight-port homodyne jnj) although we see here a 
3 x 3 - port is enough. For a general probe field we will find that triple coupler 
homodyning allows access to all possible phase space density K$p, provided the 
corresponding field pp can be generated. The use of a lower number of modes 
results in an easier and less noisy experimental implementation. 

A linear, symmetric three-port optical coupler is a generalization of the 
customary lossless symmetric beam splitter. The three input modes dj, i = 
1, 2, 3 are combined to form 3 output modes bj, j — 1, 2, 3. In analogy to lossless 



beam splitters which are described by unitary 2x2 matrices 12 , any lossless 
symmetric triple coupler is characterized by a unitary 3x3 matrix [ p"3"[ [Til of 
the form 




1 1 

exp{»f } exp{-zf } 1 , (5) 
cxp{-i4f} exp{i^L} 
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where each matrix element Tij represents the transmission amplitude from the 
i-th input port to the j-ih output port, that is bj = Ylk=i Tjk a k- 

Such devices have already been implemented in single-mode optical fiber 
technology and commercial triple coupler have been available for some time Jl5| . 
It is known |l6[ ] that any unitary M-dimensional matrix can be factorized into 
a sequence of 2-dimensional transformation. Thus any triple coupler can be 
implemented by discrete optical components using 50:50 beam splitters and 
phase shifters only |l3|, which is schematically displayed in Fig. [I] together with 
our notation for general linear three-port couplers. Experimental realizations of 
triple couplers has been reported for both cases, the passive elements case and 
the optical fiber one [|l3| [l4| . 

Let us now consider the measurement scheme of Fig. |^. The three input 
modes are mixed by a triple coupler and the resulting output modes are subse- 
quently surveyed by three identical photodetector. The measured photocurrents 
are proportional to 7„, n = 1, 2, 3 given by 

1 . 3 2n 

In = b^bn = exp{i6»„(7 - k)} a\ai , 6 n = — (n- 1) . (6) 

k,l=l 

After photodetection a Fourier transform (FT) on the photocurrents is per- 
formed 

1 3 

T s = FT(7 1 ,7 2 ,7 3 ) = - ? =^7„exp{-^ n ( S -l)} , s = l,2,3. (7) 

* n— 1 

This procedure is a straightforward generalization of the customary two-mode 
balanced homodyning technique. By means of the identity 



S 3 (s-1) = i^exp|iyn(s- 1)1 

n—1 * 



(8) 



for the periodic (modulus 3) Kronccker delta £3, we obtain our final expressions 
for the Fourier transformed photocurrents 

I\ = —j= |a}ai + a\a 2 + aja 3 | , (9) 
I2 = — J= |ai«2 + c4 a 3 + «3«l| j (10) 
and 2 3 = — |a}a 3 + a\a\ + aja 2 | . (11) 

X\ gives no relevant information as it is insensitive to the phase of the signal 
field, whereas I 2 and I3 are hermitian conjugates of each other and contain the 
relevant information in their real and imaginary part. 
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In the following let us assume a± is the signal mode and ai is fed by a 
highly excited coherent state \z) representing the local oscillator. Since the 
local oscillator serves as the reference phase (ip — 0), we have \z) = \\z\). For 
large z the output photocurrents are intense enough to be easily detected and 
can be combined to give the reduced photocurrents 

ft = v / 3^p = a 1 (o) + a 3 (o) + o[^] 

% = v / 3^^ = a 1 (-7r/2)-a 3 (-7r/2) + od], (12) 
2i\z\ \z\ 



which we refer to as the triple homodyne photocurrents. In Eq. ( p_2| ) a(ip) = 
\/2{a s e l '- p + ae~ lip ) denotes a quadrature operator of the field. 

Each experimental outcome from a fritter homodyne detector is a pair 
(1/1 , j/2) of real numbers from the joint measurements of y± and y 2 . The cor- 
responding two-dimensional probability distribution of such outcomes is given 

by 

P(jn,to)= I P I pe-^-^Z(\ u \ 2 ), (13) 

which is a double Fourier transform of the characteristic function H(Ai, A2) of 
the device, namely 

~(\ 1 ,\ 2 ) = Tr{{p 1 ®p 3 )e^+^} . (14) 

Inserting Eq. ( |l2| ) in Eq. (|l4]), the characteristic function can be rewritten in 
terms of the complex variable 7 = (A2 + i\x)/2 yielding 

2(7,7) - Tr{( / } 1 ®p 3 )e^-^ ai e-^l+^sj 

= Xo|i(7,7)Xo|3(-7,-7) • (15) 

It factorizes into two characteristic functions xo> see Eq. (Q), one for each mode. 
This and Eq. ( |l3| ) - changed to variables 7 - gives P as an integral over the 
complex plane 

P(Vi,V2)= [ ^e^^)^^ + ^) X o|i(7,7)Xo|3(-7,-7) ■ (16) 



We find that P is a convolution integral reminiscent of K$p in Eq. (g) , together 
with Eq. (0) this leads to our main result: 

P{yuyi) = #13(3/1 - W2,yi + iy%) ■ (17) 

It shows that the outcome probability distribution from a triple coupler homo- 
dyne detection gives the desired phase space density of Eqs. (||) and (||), the 
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mode 03 represents the probe mode P whereas a\ is the signal mode S. When 
the probe is left in the vacuum state, Eq. (|l5|) reads 



3(7,7) 



Trjpx® |0)(0| e 

Xo|i(7>7)e~* 1712 




X-i|i(7)7) > 



(18) 



namely, the triple coupler homodyning device becomes a Husimi Q-function 
detector. Note that customary eight-port homodyne thus mimics this scheme 
with one more input and output. 

We end the paper with Fig. |], where we illustrate some examples of phase 
space ityropensities accessible by triple coupler homodyne detection. A co- 
herent state of real amplitude (3 = 1 is chosen for the signal: Feeding vac- 
uum into the probe mode (13 the Husimi Q-function K\p\\Q\(a, a) = Q(a,a) = 
l/7rexp{— \a — f3\ 2 } is obtained, if a squeezed vacuum |0, r) is used instead, the 
distribution K\p\\n tr \(a, oc) — 1/tt| (a, | 2 results, and a number state \n) gives 
the distribution K\p)\ n )(a, a) = 1 / 'Tr\(n\& '(a)\(3}\ 2 . 
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Figure 1: Schematic diagram of a triple coupler and its realization in terms of 
50:50 beam splitters (BS) and phase shifters 'ip\ In order to obtain a symmetric 
coupler the following values has to be chosen: ipi = arccos(l/3) 
and ip2 = <Pi/2- 



Figure 2: Outline of triple coupler homodyne detectors: The hexagonal box 
symbolizes the electronically performed Fourier transform (FT). 



Figure 3: Examples of phase space densities measurable by the triple coupler 
homodyne detection. For a coherent state of amplitude (3 = 1 we find (a) 
the Husimi Q-function obtained if the probe mode az is the vacuum, (b) the 
distribution which results if the probe is a squeezed vacuum |0, r) with the mean 
photon number sinh 2 r = 1 , and (c) the distribution which results if the probe 
is a number state \n) with n = 1. 
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